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Abstract 

We study Ihara’s zeta function for graphs in the context of quivers arising from gauge 
theories, especially under Seiberg duality transformations. The distribution of poles is 
studied as we proceed along the duality tree, in light of the weak and strong graph versions 
of the Riemann Hypothesis. As a by-product, we find a refined version of Ihara’s zeta 
function to be the generating function for the generic superpotential of the gauge theory. 
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1 Introduction 


Quiver gauge theories have over the last two decades become fruitful in the cross-fertilization 
between physics and mathematics. As quantum held theories, especially those with supersym¬ 
metry, they are archetypal of those arising from string theory and phenomenology; as hnite 
graphs, they crystallize the underlying geometry and algebra, especially those with Calabi- 
Yau properties. More recently, the dialogue has extended to number theory, in particular to 
algebraic numbers and dessins d’enfants as well as to hnite helds. 

In [18], computations of the zeta-function for hnite graphs [1] was initiated for quiver 
gauge theories bearing in mind the ultimate hope of hnding the relation between properties 
such as whether the graph satishes the analogue of the Riemann Hypothesis and the algebraic 
geometry of the Calabi-Yau moduli space of vacua. Perhaps the most remarkable action on 
quivers is mutation, where in the physics literature this is a guise of Seiberg duality [19] and in 
the mathematics literature this is realized as cluster transformation [29] . It is a beautiful fact 
that these were discovered independently, one as a duality between quantum field theories 
and another as an isomorphism of quiver representations (and Calabi-Yau geometry). Thus it 
is natural to investigate the properties of the zeta function under such a key transformation. 

There is therefore a trio of conversations: the physics of the quiver gauge theory, the 
geometry of the moduli space of representation (equivalently the vacuum moduli space of 
the held theory), as well as the number theory of the zero/pole structure of the graph zeta- 
function of Ihara. It is this trio that we wish to analyse. 

We begin, in §2, with a recapitulation of the graph zeta function, for directed and undi¬ 
rected cases - though most of our quivers are directed - emphasizing on concepts such as 
regularity and Ramanujan. In parallel, in §3, we briehy summarize the rudiments of the 
quivers and associated Seiberg dualities in the gauge theory. After examining in detail case 
studies of two of the most famous quiver gauge theories, namely those arising from the world 
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volume physics of D-branes in the back-ground of affine Calabi-Yau threefolds as cones over 
and CP^ x CP^, in §4 and §5, we delve into a wealth of examples in § 6 , in order to 
see how Seiberg duality influences the Riemann and Ramanujan properties of the Ihara zeta 
function and conversely how the latter can be harnessed as a tool to explore the duality tree 
of field theories. 

In §6.1, we apply a specific graph Riemann Hypothesis for our directed quivers and 
summarize the results in Table 3. In due course, we find an interesting fact of how the 
zeta-function serves as a generating function for candidate terms in the superpotential by 
“refining” it with a multi-variable version in §7. Finally, we conclude with discussion and 
prospects in § 8 . 


2 The Ihara Zeta Function 


First, let us recall the definitions of Ihara’s graph zeta function, first defined by Ihara [1], 
studied extensively in [ 2 ,3,5-8], and introduced to the study of quiver gauge theories in [18]. 
A few graph-theoretic concepts [5], most of which are self-explanatory, are needed before we 
give the full definition for the Ihara zeta function, which we include here for completeness 
and to set our nomenclature. 


• An undirected graph G = (R, E) is a finite non-empty set of vertices V and a finite 
multiset E of undirected edges or unordered pairs of vertices. A graph is simple if there 
is no loops, f.e., no edges of the form {u,u) Vtt e R and there is only a single edge 
between any two vertices. 

• A directed graph G = {V,E) is a finite non-empty set of vertices V and a finite 
multiset E of arrows or ordered pairs of vertices. For an arrow e = {u,v), u is defined 
to be its origin o(e) and v to be its terminus t{r) (sometimes these are also called head 
and tail, respectively. However, we will reserve the word tail for a usage to be introduce 
shortly). Its inverse arrow e is defined to be (u, u) by reversing its origin and terminus. 
In general, directed graphs do not need to have both arrows and its inverse arrows both 
present in E. We make this distinction between undirected and directed graphs because 
the form of the zeta-function, as we shall see, is sensitively dependent thereupon. In 
the context of the gauge theories which we will soon study, we will primarily, for the 
sake of chirality, study directed graphs. 

• A cycle c of length n in graph G is a sequence c = (ei, 62 ,..., e^) of n arrows in G 
such that t{ei) = o(ei+i) for 1 < i < n — 1 and t{en) = t{ei). Cycle c is said to have 
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backtrack if ej+i = e* for 1 < i < n — 1. In addition, c is said to have tail if e^, = ei. 

• The in-degree (respectively out-degree) of any vertex of a graph is simply the number of 
in-coming (respectively out-going) arrows, on the other hand, the undirected degree of a 
vertex is simply the number of length -1 walks starting from this vertex using undirected 
edges only. 

• The n-multiple of a cycle c is the cycle formed by traversing c n times. A cycle is 
called primitive if it is not some n-multiple of some other cycle for n > 2. For a cycle 
c = (ei, 62 ,..., e„), the equivalence class [c] is defined to be the cyclic permutations 

[c] = {(ei, 62 , . . . , 6 n), ( 62 , 63 , ... , Cn, 61 ), . . . , (Cn, 61 , ... , 6 n-l)}, 

which simply means cycles are equivalent up to choices of initial and terminal vertices. 
Therefore a prime in a graph is primitive cycle that is non-backtracking, tailless and 
not a n-multiple cycle. 

• We are primarily concerned with graphs which are finite, i.e., finite number of nodes 
and edges) and connected, i.e. every node can be reached by traversing along some 
combinations of paths. 

• The adjacency matrix for a graph with n nodes is an n x n matrix and has its 
{i,j)-th. entry specifying the number of undirected edges from node i to j with i- 
th diagonal entry being twice the number of self-adjoining loops on i-ih. node. We 
emphasize that the adjacency matrix we use here has this convention for the diagonal; 
this has deep implications in tune with Cartan matrices, Frobenius eigenvalues and 
ADE classifications of Dynkin diagrams [14,15]. Moreover, any undirected edges in 
a graph G can be replaced with a bi-directional arrow, thereby assigning standard 
directions onto G. 


Now, we are readily for the central definition of our paper. The Ihara zeta function, also 
called graph zeta function, was initially defined [1] for graphs G that are finite, undirected, 
connected and tailless. Nevertheless, the graph is allowed to have multiple edges between 
nodes as well as loops (length 1 cycle). With these conditions, the Ihara zeta function is 
defined as: 

cgw= n (i-^'<''»)"\ (1) 

[P]ePrime Cycles 

where the infinite product is over all prime equivalence classes and ^([T’]) is the length of 
the prime cycle. This definition is clearly motivated and parallels well that of the famous 
Riemann zeta-function, whose Euler-product is C{z) = 0 (1 ~ P~^) ^ 

pEPrimes n>l 
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Indeed, as with all zeta-functions, the interplay between the expression as a product and 
as a sum is that between primality and integrality. This is also the case with the graph zeta 
function. It is shown in [8] that the Ihara zeta function takes the following closed form for 
any partially directed graphs, in which both arrows and edges are allowed [2,3]: 

(1 _ ^2)-Tr(Q-I)/2 

^ - Az + Qz^ + Pz^y 

where the following clarifications on notation are understood: 


• This is the general equation for partially directed graphs, but we only have fully directed 
graphs for particular gauge theories. 

• A stands for the full adjacency matrix with its entry aij specifying length 1 walk from 
i-th node to j-th node using either an edge or an arrow. 

• P is the directed adjacency matrix whose entries are composed entirely of arrows. 

• Q is the matrix for undirected degrees. The i-th undirected degree is the number of 
length 1 walks from i-th node to its neighbours using undirected edges only. Specifically, 
it turns out that the exponent in the numerator stands for the number of vertices minus 
the number of edges, or simply 1 — r with r being the rank of fundamental group of the 
graph. For nodes with self-adjoining loops, the undirected degree is counted as 2. 


Lastly, the relation between Euler-product and determinant expression of of Ihara zeta func¬ 
tion is discussed in more detail in Appendix A. 


2.1 Poles, Regularity and Ramanujan 

Some of the initial motivations for defining the zeta function of graphs are, of course, number 
theoretical, and in particular in issues such as finding analogues of the Riemann Hypothesis 
(cf. [16] for a recent equivalent restatement of the Riemann Hypothesis in string theory). 
Here, let us highlight some salient features of Ihara’s zeta function. 

Now, a graph is called a {q + l)-regular graph if all of its nodes are connected to other q + 1 
nodes through length one walks. In addition, a directed regular graph must satisfy a stronger 
condition that all the nodes must have their in-degrees and out-degrees equal (see page 29 
of [17]). A (gr-|- 1)— regular graph is called Ramanujan if the maximum of the absolute value 
of eigenvalues of the adjacency matrix A, excluding q + 1 itself, is bounded by 2y^. In short, 

maxjjAj : A G Spec(A), JA] / g -|- 1} < 2y/q . (3) 
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While the original motivation to study Ramanujan graphs was because they have maximal 
gaps in their spectrum, it turns out, as we now see, that they play a key role in a graph version 
of the Riemann Hypothesis. The Ihara zeta function of a (g + 1)—regular undirected graph is 
said to satisfy Riemann Hypothesis if, in complete analogy with the number-theoretic case, 
for 0 < Re(s) < 1, the zeros of Cg{q~^)~^ He on the Re(s) = ^ line. Here, as is customary 
with definitions of zeta-functions, we define the exponentiated variable 

.2 := q-^ . (4) 

Note, however, the parameter q is natural for regular graphs whereas for irregular ones more 
work is needed to extract a similar quantity. In addition to the above definition for Riemann 
Hypothesis for undirected regular graphs, it also can be shown [5] that a (q + 1)-regular 
undirected graph satisfies Riemann Hypothesis if and only if it is Ramanujan. 

Apart from the aforementioned property, Ramanujan graphs have diverse connections 
with various fields, such as expander graphs in Communication Network Theory that revolves 
around extremal problems. Number Theory, Representation Theory and Algebraic Geometry. 
For more comprehensive surveys, we refer the readers to [5,9,10,12,13]. A more recent survey 
concentrating on connection between Graph Theory and automorphic representation can be 
found in [11]. 

In more generality, we also have a definition of the Graph Theory Riemann Hypothesis 
for irregular undirected graphs as follows: 


1. Cg{z) is pole free for 

Rg < \z\ < '/Rg, ( 5 ) 

where Rq is the radius of convergence of Cg{z)^ *-e., the position of nearest pole to 
origin. 


2. There is a weaker version of the above, requiring Cg{z) to be pole free in 


Rg < k! < 


1 

7 ^’ 


where q is the largest degree of vertices. 


( 6 ) 


These two definitions are motivated by the fact that if z is substituted by Rq, all poles 
of Cciz) are then located within the “critical strip”, 0 < Re(s) < 1, with poles at s = 0 (i.e., 
z = 1) and s = 1 (i.e., z = Rg)- Importantly, as is central to the study of any zeta function, 
there is an underlying functional equation. It can be shown [5] that the functional equation 
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for Ihara’s zeta function for undirected graphs can take the form 

Kg{z) := (1 - z2)-i+i(i _ q^z^)"HG{z) = (-1)"Ag(-), (7) 

qz 

where r is the rank of fundamental group with n being the number of vertices. With substi¬ 
tution z = q~^, it is obvious that the above equation is symmetric with respect to s = ^. 

Since irregular graphs do not have a functional equation relating f{s) and /(I — s), it is 
natural to define a pole free region for ^ < Re(s) < 1, therefore motivating (5). The weak 
version in this sense simply shrinks the size of the pole free region. 

It should be noted that the formulation of Riemann Hypothesis in terms of Ihara zeta 
function is based on the fact that the adjacency matrix of an undirected regular graph is 
symmetric and we can find a functional equation f{s) which is symmetric with respect to 
s = ^. However, for the case we will consider shortly, especially in the context of Seiberg 
duality in A7 = 1 supersymmetric gauge theories, we will be primarily concerned only with 
directed graphs. Since the adjacency matrix in these theories are rarely symmetric and the 
form of Ihara zeta function is very different from the undirected case, it is not possible to 
construct a functional equation similar to that of (7), thus we can not construct a direct 
analogue of Riemann Hypothesis for directed graphs. 


3 Four-Dimensional Quiver Gauge Theories 


Having recapitulated the relevant ingredients from the study of graphs and Ihara’s zeta 
function, we now turn to our protagonist, the finite (generically directed) graphs which are 
realized as quivers associated to four-dimensional, N = 1 supersymmetric gauge theories. 
Briefly, quivers are finite, labelled, directed graphs. In addition to these properties, they are 
also allowed to have loops, bi-directional arrows and self-adjoining loops. On the vertices we 
have gauge groups and the corresponding gauge theory has product gauge group in the from 
of n77(fVj) with f-th vertex contributing a single gauge group factor. In the infra-red, the 
U{\) subgroups of each of the U{Ni) factor become frozen and we have an effective gauge 
group of products of special unitary groups. However, from the quiver representation point 
of view, it is important to retain the U{Ni) as labels of the nodes. The matter content of 
the theories are represented as bi-fundamentals between gauge groups, so each arrow from 
f-th vertex to j-th vertex transforms as {Ni,Nj) representation of the factor gauge group 
SU{Ni) X SU{Nj). Gauge theories with only undirected edges in their quivers are non- 
chiral and i-th vertex should be associated to a vector space with arrows being maps in 
Hom(C^SC^^). 
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Furthermore, the matter content of the theories mnst be anomaly free. This condition 
ensures that the corresponding quantnm field theory is well-defined. When translated into 
quivers, this constraint has the following form: 

{oij - aji)ni = 0 , ( 8 ) 

where Uij is the full adjacency matrix of the quiver with n* being the vertex rank vector 
encoding the list of integers W- This condition will become important when we discuss the 
zeta function of our quivers, and certain Diophantine equations dictating their possible rank 
assignments. 


3.1 Seiberg Duality and Cluster Mutations 

Modern studies of AA = 1 supersymmetric gauge theory lead to Seiberg’s discovery of IR 
duality between two QCD-like theories [19], where the dnals have same Nf fundamental 
chiral flavors of quarks, but different gauge groups SU{Nc) and SU{Nf — Nc). The duality 
predicts that the two theories have the same supersymmetric moduli space and ’t Hooft 
anomaly matching condition. In extending the duality to other M = 1 theories [20-22] , the 
most suggestive method is via the Hannany-Witten suspended brane construction [23]. In 
snch constructions, we obtain a D-dimensional gange theory by extending a set of Dirichlet 
branes in D dimensions, i.e. the embedding in the transverse 10 — D dimensions determines 
the spectrum and other properties of the theory. Therein, the D-branes are strings probing the 
transverse space, with one or both ends on the NS5-branes. If we take Nc finite length strings 
between 5-branes, a U(Nc) pure world-volume gauge theory is produced. By introducing Nf 
strings with semi-infinite length, Nf quarks are produced. 

In this stringy context, Seiberg duality is obtained through moving one 5-brane to ex¬ 
change its position with others, while reversing the orientation of the finite strings. Within 
the context of Nf > Nc, the movement and reconnection of strings prodnce a dnal theory with 
Nf semi-infinite and Nf — Nc strings [24]. Moreover, since the inverse conpling is identified 
with the length of the string, the duality arises if we vary the gauge coupling through infinity 
via variation of a moduli field. Since the embedding of gauge theories into string theory is 
not unique, the field-theoretic duality is interpreted as different brane configurations that 
give dual low-energy physics. 

Before we discuss Seiberg duals for specific theories, it is worthwhile to look at the rules 
for Seiberg transformation on gauge theories. We summarize that the prototype [19] is a pair, 
within the conformal window < Nf < 3Nc, (1) direct electric theory with Nc colours 


with Nf flavours and (2) dual magnetic theory with Nf 


Nc colours also with Nf flavours: 



SU{Nc) 

SU{Nf)L 

SU{Nf)ji 

Q 

□ 

□ 

1 

Q' 

□ 

1 

□ 


W = 0 



SU{Nf - Nc) 

SU{Nf)L 

SU{Nf)R 

q 

□ 

□ 

1 

q' 

□ 

1 

□ 

M 

1 

□ 

□ 


W = 

= Mqq' 



(9) 


In the above, the quarks Q and Q' are transformed to the dual quarks q and q' and a Seiberg 
dual meson M together with superpotential Mqq' is generated. It was recognized in [25-28] 
that, when applied to one node of any M = 1 quiver, Seiberg duality is the following graphical 
rule: 


Nc 


Nf-Nc 



A truly remarkable fact is that around about the same time, mathematicians have indepen¬ 
dently noticed the importance of this graphical rule in the study of cluster variables [29] and 
the above was called cluster mutation (though at the time, the role of the superpotential was 
not yet appreciated in the mathematics). That Seiberg duality as an equivalence of quan¬ 
tum theory should be the same as cluster mutation generating an equivalence in the derived 
category of coherent sheafs on Calabi-Yau manifolds is a deep result in mathematical physics. 


From an algorithmic point of view, one can think of the duality/mutation as the following 
transformation, for a quiver gauge theory with adjacency matrix aij and n vertices or gauge 
group factors, on a^j and rank vector n*. Indeed, there are n choices to perform Seiberg 
transformation at any stage and the arrows and ranks will be changed while n will always 
remain fixed. Suppose, we dualized on node io, then the transformation is: 


1. Define Iin to be the nodes having arrows coming into io, lout to be those having arrows 
coming from io and I no to be those unconnected to io. 

2. Change the rank of node io from Nc to Nf — N^ where Nf = ®*.*o-^* ~ 

J2ieiout This changes the rank of the gauge group of the dualized node. 

3. = Qji if either i,j = io. In field theory context, this translates to the statement 
that the quarks of SU{Nf — Nc) gauge group are in complex conjugate representation 
to the quarks of original SU{Nc) group. 

4. — aifjaajSio for ck G lout and f3 G lin- This is equivalent to adding Seiberg 
meson into the dual theory. Any bi-directional arrow corresponds to a quadratic mass 
term and can be integrated out. 
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4 Illustrative Example: Cone over 


Having presented the necessary mathematics in terms of graph zeta functions and associated 
statements pertaining to their Riemann Hypotheses, as well as the physics in terms of four¬ 
dimensional M = 1 quiver gauge theories and their Seiberg duality, it is illustrative to start 
with an archetypical example to see the interplay of the two sides. 


4.1 dP^: Cone over 


Perhaps the most well-studied quiver is the one presented in Figure 1. If we write the 3 
multi-arrows as Xi,Yi,Zi, it has an accompanying superpotential W = . This 

quiver gauge theory comes from taking the quotient of AA = 4 super-Yang-Mills theory by a 
Za-subgroup of SU{3) to preserve M = 1 supersymmetry. Geometrically, since the moduli 
space of the parent J\f = A SYM is simply the affine Calabi-Yau 3-fold C^, the moduli space 
here is jlLj, as a complex cone over (and whose resolution can be seen as the total space 
of the anti-canonical bundle over P^). Thus, historically, the theory is called the P^ or 
dPo theory. 


N} 



iV3 


Figure 1: Quiver for P^, vertices have SU{Nc) factor gauge groups and the number of arrows 
between each pair of vertices is 3. 


4.2 Adjacency Spectra and Poles 


Since this quiver is fully directed, its Ihara zeta function has a particular simple form [3,18]: 


Cg{z) 


1 

Det(/ — Az) ’ 


( 10 ) 
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where A is the adjacency matrix for digraph G. We see therefore that the reciprocal of 
zeta function for directed graphs is none other than the characteristic polynomial of the 
adjacency matrix A, 

XaW = Det(A/ - A) ~ Det(/ - . (11) 

Therefore the behaviour of poles (those z that make Cdz) singular) in Ihara zeta function is 
inversely proportional to the distribution of eigenvalues (those A that make xa(A) vanish) of 
adjacency matrix for each Seiberg dual quiver. 


Let us call the quiver in Figure 1 Qo and hence the adjacency matrix and the Ihara zeta 
function are 

/ 0 3 0 


= 


0 0 3 
V 3 0 0 


CQo(^) = 


1 


1 - 27^3 


( 12 ) 


The above will be our “ 6 asic case” from which we will perform repeated Seiberg duals on 
different nodes. Using rules from Section 3.1 to dualize on node N^, we have following results 
as presented in Figure 2. 


N} 


Nj - iV 1 




iV3 


Figure 2 : Quivers for and its Seiberg dual on node N^. 


The reciprocals of their zeta functions are readily seen to be 1 — 27^3 and 1 — 542;3 
respectively. Hence the poles lie on the lines of third roots of unity with \z\ = 1/-^|c 3 | where 
C 3 is the coefficient of term in the reciprocal. Let us perform duality to a few levels, picking 
any of the three nodes each time. This has been well-known [32,33,35] to exhibit a dendritic 
behaviour and can be drawn as a duality tree where each node in the branching corresponds 
to a new quiver (i.e., this is a tree in the space of theories). For reference, the tree for dPo 
is included in part (a) of Figure 3, where the same colours correspond to the same quiver, 
reflecting the 3-fold symmetric in the problem. 
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Figure 3: duality tree and the position of the poles on the complex plane for the lhara zeta 

function for all the dual quivers, up to 5 levels. NB: Colors of leaves in part (a) and the colors of 
dots in part (b) are not related. As we delve deeper into the duality tree, the lengths of poles, on 
a whole, get shorter, as the poles accumulate toward the origin. This can be seen from Table 1. 


Now, let us compute the poles of the zeta function, for the various dual phases. Numeri¬ 
cally, we see that their length (i.e., distance to the origin on the complex plane) are as shown 
in Table 1. In fact, all the Seiberg duals have their zeta function poles sitting on the lines 


Level in duality tree 

All possible pole lengths 

1 

1/3 

2 

0.26 

3 

0.15 

4 

0.083, 0.050 

5 

0.044, 0.016, 0.014, 0.0084 

6 

0.023, 0.0048, 0.0040, 0.0024, 0.0014, 0.0012, 0.00080, 0.00043 


Table 1: Lengths of poles in first 6 levels. By length we simply mean the distance of the pole to 
the origin in the complex plane. 

of cube roots of unity for quiver (see Figure 3) and this can be derived analytically as 
follows. 

If we denote (tt, u, w) as the number of bi-fundamental arrows between the 3 nodes, a tree 
of branching integer triplets can be obtained by successive dualization. These triplets in turn 
describe all possible solutions of the P^ quiver under Seiberg duality. Furthermore, for all 
dual solutions to be anomaly free, the Diophantine equation [30-35]: 

+ v‘^ + w‘^ = uvw (13) 
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has to be satisfied. This is the famous Markov equation and it is well-known that all its 
solutions can be generated from the basic solution (3,3,3) by performing ad infinitum the 
following three transformations [37] (the proofs are gathered in Appendix B for the readers’ 
reference), 

{u,v,w) —>■ (u,v,uv — w), 

(u,v,w) —>■ (u,uw — v,w} , 

(u,v,w) —)■ (vw — u,v,w). (14) 


One can see that the above transformations are similar to Seiberg transformations per¬ 
formed on quiver arrows (see rule 3 and 4 at the end of Section 3.1). To be specific, the 
Seiberg transformation rule on arrow numbers are, 

{u,v,w) —)• {—u,—v,w — uv) , 

{u,v,w) —)• {—u,v — uw,—w) , 

{u,v,w) —)• {u — vw, —V, —w) , (15) 

which differ from (14) only by a common factor —1. This common factor can be interpreted 
as a reversal in the arrow directions. So we can say, a Seiberg transformation generates a 
new quiver whose arrow number combination is another solution to Markov equation which 
means this new quiver theory is automatically anomaly free. Since (14) generates all possible 
solutions, we can conclude that Seiberg transformation generates all possible anomaly free 
dual quiver theories. 


On the other hand, since Seiberg transformation reverses arrows, if we choose the basic 
quiver to be Qq whose cycles are all anti-clockwise as shown in Figure 1, then all its descen¬ 
dants have their cycles being either clockwise or anti-clockwise. Due to this fact all adjacency 
matrix A will have one of the following two forms. 


/ 0 

u 


f ° 

0 

u \ 

0 

0 

V , 

v 

0 


\ w 

0 

oj 

[ 0 

w 

0.) 


In either case we have the same Ihara Zeta function C{^)^ 

(■”^( 2 :) = det(/ — Az) = 1 — {uvw)z^. 

Thus the poles of Zeta function are 

1 1 

^uvw ’ \fumu 


(16) 


(17) 


(18) 
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This explicitly demonstrates that all poles after Seiberg transformation lie on lines of third 
roots of unity (see Figure 3). Moreover, since we have all the information on poles for all 
Seiberg transformed quivers in this particular case, according to equation (11), we can see 
that all adjacency spectrum of each transformed quiver should also lie on the same lines as 
the poles. 

It is interesting to remark that since the original quiver already has arrows between 
every pairs of nodes, Seiberg transformation does not connect any previously unconnected 
nodes to generate any new loops of length other than 3. As a consequence of this, the 
reciprocal of Ihara Zeta functions of all duals only have cubic terms in the complex variable. 
Indeed this property of Ihara Zeta function is rather general and how it gathers information 
in terms of simple cycles in a graph will be discussed in detail in Section 7. 


4.3 Ramanujan Condition for and its Seiberg Duals 

Due to the fact that Ihara zeta function only has the analogue of Riemann Hypothesis for 
undirected graphs, i.e. for {q + l)-regular undirected graph G, Cci^) satisfies the Riemann 
Hypothesis if and only if G is Ramanujan, and quivers under consideration being fully directed 
graph, it is only possible to check if the graph is Ramanujan. It should be noted that 
Ramanujan condition for directed graph has an extra constraint. In addition to the definition 
in section 2, the adjacency matrix of directed graphs should be diagonalizable by unitary 
matrices. See [36] for more details in the definition. 

First of all, a quiver being Ramanujan requires this quiver to be a regular graph (that 
is, all vertices have the same in-degree and out-degree d). Since only the singular solutions 
to (13) can have same nnmbers and only the basic solution has its three numbers the same, 
obviously only the basic quiver is a regular graph. Thus none of the other quivers can be 
Ramanujan. We then check whether this basic qniver, denoted earlier as Qo; is so. We see 
that here all the eigenvalnes A have |A| =3, which means Qq cannot be Ramanujan and thus 
violates the graph Riemann hypothesis. 


5 Example: Cone over the Zeroth Hirzebruch Surface 


Another well-stndied example in both the physics and mathematics literature (independently) 
is the quiver which corresponds to the cone Fq over the zeroth Hirzebruch surface, Fq := 
X F^, which is a toric variety. The archtypal Fq quiver (we shall call it the basic quiver) 
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can be visualized in Figure 4. 


0 


2 



Figure 4: Basic quiver for x vertices are labeled with 1,2,3,4. 


A series of quivers - we shall call them general Fq quivers - are generated by Seiberg 
dualizing on different vertices for a number of times and they are pictured in Figure 5; 
indeed, this is the most general form of a 4-node quiver. We use the convention that a 
positive integer denotes a multiple of arrows while a negative means the same, but in the 
reverse direction. Perhaps the most famous dual quiver to the (a, b, c, d, e, /) = (2, 2, 2,0,0, 2) 
basic case is {a,b,c,d,e,f) = (2,—2,—2,4, 0, 2) on dualizing on node 1, because both these 
quivers afford rank vectors which are all 1, i.e., they are both toric. 



Figure 5: General quiver emerging from dualizing the basic quiver. 

The antisymmetric adjacent matrix q can be read off from this general quiver as 

^0 c —e 

—c 0 o 

n 

e —a (J 

\ b d -f 

On the one hand, to guarantee that the kernel of q is not null so that one could have a 


-b ^ 
-d 
f 
0 
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sensible rank-vector according to anomaly cancellation in (8) we should impose the following 
condition (cf. [35]) 

det q = [ab + de — c/)^ = 0. (20) 

On the other hand, as a consequence of gauge anomaly cancellation, a Diophantine equation 
is found [34], 

+ b^ + (? + + bed + bef = ace + adf -|- abef. (21) 

Physically one would expect every Seiberg dual theory of the basic quiver is an anomaly free 
theory, i.e., each edge number combination in these dual quivers forms a solution to the above 
two equations (20) and (21). 

Yet, we seem to find that the converse proposition also holds. That is, all positive even 
solutions to these two equations can be generated by Seiberg duality starting from the basic 
quiver. By “positive even solution” we mean all variables in this solution take positive 
and even integer values. Indeed this a very complex set of multi-variable and non-linear 
Diophantine equations and we need to use (20) to get rid of one of the variables. Firstly, it 
is easy to see that no solution can have c = 0. Suppose we have a solution with c = 0, then 
a or 6 and d or e would have to be 0, the only possible non-trivial combination has to satisfy 

6 = c = e = 0, -b -I- /^ = adf , (22) 

which reduces to the quivers discussed above and is not the case we want to discuss here. 
With this restriction, we have / = {ab + de)/c from (20) and equation (21) then reduces to 
a 5-variable Diophantine equation, 

n 2 l2 I 2 2 I f 2 2 i 2 t2 i 2 2 i y2 2 i 4 2 f 2 2 i l 3 y 

0 = ab+ac+bc+cd-\-ce+de+c—abc+ be d 

—a(?e -b ab^ce — a^bed + 2abde — acd^e -b bede^ — ab(?de. (23) 

Instead of solving this equation analytically which is beyond our reach we seek to hnd as 
many solutions as possible with computer power and use these solutions to test our conjecture. 
The basic algorithm is as follows. 

1 Generate all combinations of (a, b, c, d, e) with any of the variables being a positive even 
integer within the range [0, Ai]; 

2 Test each combination and select those that solve equation (23). 


The result is, we have found all solutions with N = 800 and all of them can be found in 
the theories of general quivers. This is the analogue of the statement that Seiberg duality 
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generates all solutions to the Markov equation and it would be nice to have an analytic proof 
of this fact. For more detailed reference of our argument, we have put all results and the 
source code at GitHub^. 


5.1 Ramanujan Condition for Fq and its Seiberg Duals 


As with the dPo case, we can check explicitly if all the dual quivers generated from the basic 
X quiver are Ramanujan. First, let us check the basic case. The adjacency matrix there 
is 


^ 0 

2 

0 

0 ^ 

0 

0 

2 

0 

0 

0 

0 

2 


0 

0 

0 j 


(24) 


and the poles are computed explicitly to be ± ^ and ^ 2 ^ and the modulus of all eigenvalues 

are thus seen to be |A| = 2. Therefore the basic case is again not Ramanujan. 


For the Seiberg duals generated from P^ x P^, it is also checked in 23 that the only regular 
graphs are those with the same Ihara zeta function, (1 — This is equivalent to say 

that the only regular graphs in P^ x P^ and its duals are P^ x P^ itself and the one that has 
all its arrow in opposite direction compared to P^ x P^. 


6 Zeta Poles and Gauge Theories: A Plethora of Examples 


With our experience now in the dPo and Fq examples, we can proceed to study the plethora 
of quiver gauge theories known to the literature. We will present pole plots generated from 
Seiberg duals of Fq, dPi, dP 2 and dP^ quivers. Since an analytical result for poles of dPo 
and its duals was obtained in section 4.2, it is omitted here. These Hirzebruch and del 
Pezzo theories [25] have since become canonical examples of studying (toric) quiver gauge 
theories and for the reader’s convenience we tabulate the relevant (basic toric) quivers and 
their associated Ihara zeta functions here (cf. [18]) in Table 2. 


^https://github.com/dayzhou/Seiberg-Diophantine. 
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While in the above table we have distinguished the toric phases because they are most 
popular, there is no reason to restrict to them here. Indeed, we compute all the zeta functions 
along the duality tree for all the above geometries and study the pole structure thereof at 
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(a) Pole plot for x P^ quiver and its duals. 


(b) Pole plot for dPi quiver and its duals. 


Im[z] Im[z] 



(c) Pole plot for dP2 quiver and its duals. 


(d) Pole plot for dP^ quiver and its duals. 


Figure 6: lhara zeta function pole plots for Fq, dPi, dP 2 and dP^ quivers and their Seiberg duals. 


once. These are presented in Figure 6. The figure is generated by dualizing on every vertex 
of the original quivers and then repeating the process for every dual quivers generated from 
the previous stage, following along the complicated duality trees (See [33]). Furthermore, 
due to the fact that the coefficients of the inverse lhara zeta functions of dual quivers get 
very large very quickly {e.g. for Don Zagier conjectured that the nth Markov number is 
asymptotically given by = igC'Vii+C’li) with C = 2.3523 ...) we only look at dual quivers 
with arrow number smaller than 80. 

The plots of Fq and dPi are very similar and it can be explained due using results from 
section 7. Since the coefficients of the term in the inverse of lhara zeta function expansion 
counts (we will return to discuss this point in detail in the following section) the number of 
length n simple cycle if there does not exist smaller disjoint simple cycles with length ni and 
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n 2 such that n = ni + n 2 - Hence in the Fq and dPi dual quivers, the highest order term in the 
inverse of zeta function is as there are only 4 vertices in total and it not possible to form 
length-5 simple cycle or higher length simple cycle using disjoint lower length simple cycles. 
Therefore the inverse of their Ihara zeta functions should all have the form 1 -|- az^ + bz‘^, 
where a,b £ Z. The solutions to these equations therefore reflects the shape of pole plots of 
Fq and dPi. Similarly, dP 2 can not have term of order higher than 5 in their expansions and 
dPs can not have term of order higher than 6 as well. 

Since all of these toric phases can be obtained from Seiberg transformation, their poles 
are also plotted in figure 6. It is also interesting to note that the second Phase of dP^ is 
the only regular graph in all of Seiberg duals of dP^. However, its adjacency matrix is not 
diagonalizable by a unitary matrix, and it is therefore not Ramanujan. 


6.1 A Numerical Experiment on Graph Riemann Hypothesis 

As we saw in the introductory sections that the Ramanujan property is intimately linked to 
the graph Riemann Hypothesis. In this subsection, we shall explore certain possibilities in 
applying both strong and weak Graph Riemann Hypothesis (RH) to our toric phases of del 
Pezzo quivers. Even though there is no immediate definition of RH pertaining to directed 
graphs due to our lack of an functional equation with reflection symmetry about s = ^, it is 
still interesting to see how certain setups of RH for undirected graphs select special theories 
from Seiberg duals of del Pezzo quivers. In this perspective, we have the following setup: 

1. The degree g -|- 1 is defined to be the maximum degree among all in-degrees and out- 
degrees of all vertices. This definition is in analogous spirit to that of Eq. (6), where 
the definition of weak RH is given for undirected irregular graph. 

2. Dual quivers are selected if they satisfy either strong or weak RH. 

3. We then look at distinct prime factors for each term in for each dual quivers 
satisfying either two versions of RH. Especially if they form any interesting prime 
sequences. 

4. Finally we look at how RH selects dual quivers in the duality tree. 


First let us consider case, since it was shown in section 4 that all of its Seiberg duals 
have poles lying on lines of cubic root of unity or having same modulus, we immediately see 
that both versions of RH are trivially satisfied. 
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Base Quiver 

Level 

RH Duals 

Prime Sequence 

pi ^ pi 

5 

1 - 96^3,1 - 1056^3, 

1 - 480^3,1 - 14432^3, 

1 - 27202:3,1 - 28320^3, 

1 - 686402:3,1 - 200736^3 

23; 2,3,5,11,13,17,41,59 

dPl 

5 

None 

N/A 

dP2 

Toric Phase 1 

5 

1 - 35^3 - 80z^ - 25^3, 

1 - 528^3 - 3082z^ - 2480^®, 

1 - 366^3 - 1594^4 - 700^3, 

1 - 984^3 - 8536z^ - 12040^3 

23: 2,3,5,7,11,41,61 

2^; 2,5,11,23,67,97,797 
23: 2,5,7,31,43 

dP2 

Toric Phase 2 

5 

1 - 35^3 - SOz^ - 25^3, 

1 - 366^3 - 1594^4 - 700^3, 

1 - 528^3 - 3082z^ - 2480^3, 

1 - 732^3 - 5328^4 - 5994^® 

23: 2,3,5,7,11,61 

2^: 2,3,5,23,37,67,797 

23: 2,3,5,7,31,37 

dPs 

Toric Phase 1 

4 

1 - 36^3,1 - 25^3 - 61^4 - 30^3 - 4z«, 

1 - 8^3 - 64^4,1 - 72^3, 

1 - 180^3,1 - 12^3 - 92z^ - 40^3, 

1 - 173^3 - 523z^ - 96^3,1 - 792^3, 

1 - 205^3 - 806z^ - 531^3 - 90^®, 

1 - 77^3 - I69z^ - 90^3 

23: 2,3,5,11,41,173 

2^: 2,13,23,31,61,523 

23: 2,3,5,59 

2®: 2,3,5 

dPs 

Toric Phase 2 

5 

1 - 36^3,1 - 25^3 - 61z^ - 30^3 - 
1 - 8^3 - 64^4,1 - 72z^, 

1 - 180^3,1 - 12^3 - 92z^ - 40^3, 

1 - 205^3 - 806z^ - 531^3 _ goz^, 

1 - 1732;3 - 523z^ - 96^3,1 - 7922;3 

23: 2,3,5,11,143,173 

2^: 2,13,23,31,61,523 

23: 2,3,5,59 

2®: 2,3,5 

dPs 

Toric Phase 3 

4 

1 - 362;3,1 - 180^3, 

1-25^3-612^-30^3-4^3, 

1 - 8^3 - 642^, 1 - 7223, 

1 - 32023 - 16882^ - 161223 - 2802®, 

1 - 1223 - 922^ - 4023, 

1 - 17323 - 5232^ - 9623, 

1 - 252023,1 - 79223 

23: 2,3,5,7,11,173 

2^: 2,23,61,211,523 

23: 2,3,5,13,31 

2®: 2,5,7 

dPs 

Toric Phase 4 

3 

1 -3023,1 - 18623, 

1 - 923 - 28^4^ 1 _ 0,j3 _ 

1 - 5123, 1 - 423 - 142^ - 623, 

1 - 623 - g^4 _ i2z^^ 

1 - 270323,1 - 1623 - 442^ - 1223, 

1 - 2223,1 - 3123 - 802^ - 3623, 

1 - 54323,1 - 4223 - 1302^ - 2823, 

1 - 523 - 9^4 _ 7^5 _ 2^6 

23: 2,3,5,7,11, 
17,31,53,181 

2^: 2,3,5,7,11,13 

23: 2,3,5,7 

2®: 2 


Table 3: Prime sequences in this table are produced for each term in C ^ expansion. Level 
specifies the level in duality tree up to which the dualization is performed. 
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We move on to the higher del Pezzo quivers and tabulate the relevant experimental results. 
In Table 3, we see that toric phases have similar prime factor sequences for each term in 
This should be expected as toric phases can be reached by Seiberg transformation from each 
other. The difference in this case can be due to the fact that the dualization started from 
different quiver and if this process is continued long enough, sequences of different toric phases 
should converge to the same one. Specifically, the dual x quivers that satisfy RH are 
those with terms only in its inverse Ihara zeta functions. This is also expected as those 
with only terms reduce to P^ case, which trivially satisfies both RHs. 

It is more illustrative to look at how RH selects quivers in the P^ x P^ duality tree, which 
we include from [33] for reader’s convenience. 



1 

•- 


© 


22 


22 
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1 

5 


2 

10 



1 


128 


65 




(a) Duality tree for x theory. (b) Quiver Rules for P^ x P^ duality tree. 

Figure 7: Duality tree and corresponding quiver rules for Fq. 

Since we have all our P^ x P^ dual quivers satisfying RH to contain terms in z^ only in 
table 3. It can be seen from figure 7 that in our setup, RH only chooses the node without 
repetition of colors. In our plot, this corresponds to a path that takes the sequence of colors 
such as blue —)> cyan —>■ purple —>■ green —>■ grey. 
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7 Ihara Zeta Function, Quivers and Superpotentials 


In this section we focus on a proof of a proposition presented below, which reinterprets 
coefficients of the inverse of Ihara zeta function in terms of simple cycles. In terms of gauge 
theories, this translates to generic superpotentials that can be generated from certain quivers. 

Proposition 7.1. For a fully directed quiver G with no self-adjoint loops, the reciprocal of 
its Ihara zeta funetion is the generator for simple loops in the sense that 

= E E (25) 

,jk}J2ik=k “ / 

where [a] is the number of ways to walk through a simple cycle in a particular vertex sequence. 


Proof. As we saw in the previous sections, the reciprocal of the Ihara zeta function for a fully 
directed graph G with no self-adjoint loops is simply 

Cq^{z) = Bet{In - Az), (26) 

where A is the adjacency matrix with entries denoted by Uij. Let us define Z as the matrix 
in the above determinant expression, 

1 - 0122 ; • • • -ainZ 

- 0212 ; 1 • • • -a2nZ 

UnlZ U.yi2Z * * * 1 

Now, Eq. (26) can be expanded as 

C ^ = DetZ = eq...i„Ziq • • • (28) 

'"'in 

where we have omitted the subscript G and argument 2 : on the left for convenience. Obviously 
the right hand side is a finite polynomial of 2 ; in which no term has its power higher than n. 
All terms have several I’s and several —aijz terms, therefore a general term with (n — k) I’s 
would look like 

~ (~^) ■ (^9) 

Written in this way, we should note that we do not sum over all i-indices any more. The 
reason is that once we have chosen some combination of {ji, ■ ■ ■ , jfc}, say {1, • • • , /c}, then 
the I’s automatically take other indices, i.e. k + 1 lo n. Since all I’s are coming from the 
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diagonal entries, which means = • • • = Z^n = 1 or m for m > fc. For this 

specihc combination, this term simplifies to 

( 1) ( f) •••**.®lii ‘ ‘ ‘ j 

which indicates only ii to are summed indices. We can then rewrite ( 28 ) as 


n 

k=Q 

Ct = (-1)‘ 

^-hi-hfc- “hhi ■ ■ ■ “ifchfc ’ 


{ii.-- ,jk} {'• 

Hi >--- ihfcl 


= (-D* E 


( 31 ) 

{h,-- dfcl {hr-- dk} 



where the first summation in Cfc sums over all possible unordered combinations of {ju •'' > Jfc} 
and the second sums over all permutations of {ji, • • • ,jk}, and it is understood that = 1 

if ii < • • • < ifc. One can easily see that cq = 1 from the above expression. Next we want to 
find out the graphical meaning of other coefficients. 

For clarity, we define some terms to make later analysis easier. 

• Simple Loop: we call a monomial a simple loop if the indices of its components can 
be written in a cyclic pattern: aijUji • • • apqUqi and each index appears exactly twice, 
once as a row index and once as a column index. We say a monomial with its indices 
written in such cyclic form is in its standard form. We use square brackets to denote 
the standard form: 

• (^ij ■ ■ ■ (kii, 

where the subscripts keep track of the vertex sequence in this loop. 

• Equivalent Loops: we say two loops are equivalent if their index sequences are different 
rotations of the same sequence given that both terms are written in standard forms. 
E.g., 012023031 = 023031012. Obviously equivalent loops have equal values. 

• Disjoint Loops: we say two loops are disjoint if their corresponding monomials use 
different (disjoint) sets of indices, e.g., 012023031 and 045054. In terms of graph theory, 
they do not share same vertices. 

Now we claim that every term in the summation of ( 31 ) can be split into several (or 1 ) 
disjoint loops. This is due to the fact that for any term in the expansion of a determinant 
each index should appear exactly twice. Then we can start from an arbitrary Ojj and select 
Ujk to be the next entry. By repeating this procedure we would complete searching for all 
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simple loops in this term with no entries left. With such statement, we can rewrite (31) as 
summation over products of disjoint simple loops, 

Cfe = (-1)^ (32) 

{iir- Jk} 

where the second summation sums over all possible inequivalent simple loop combinations 
whose lengths add up to k and whose indices/vertices ii are chosen only from the set 
{ji)■■ ■ )Jfc}• 

Now let us focus on the second summation to see if we can extract more graphical infor¬ 
mation therein. First of all we shall assume without loss of generality that the k vertices have 
already been chosen to be {1, 2, • • • ,k}. Second we would rewrite the Levi-Civita symbol e 
as products of shorter e’s each of which is associated with one simple loop. This can be most 
easily seen by a concrete example, 

£3452iai3a24«35a42a51 = £35ie42(ai3a35fl5l)(a24«42) = (^SSl [a] 135l)(^42[«]242)• (33) 

Generally, we have 

ei...iaii ■■■aki = (ei...[a]ii...i) • • • (e^... [a]...y..). (34) 

Then we can rewrite the second summation in (32) as 

T,ik=k i 

where k is the length of the i-th simple loop in one combination of loops whose lengths add 
up to k. 

We want to prove that in the above product: (I) (—l)^®ej.... = —1, (II) [a]...j-... is the 
number of ways travelling through these vertices in this simple loop. The second proposition 
can be readily understood by the fact that aij is the number of arrows from i to j which 
means there are aij ways to walk from i to j, so the number of ways going through all vertices 
are their product. The proof of the first proposition is slightly tricky and we need to take 
advantage of symmetric groups. 


Suppose we have I vertices labeled with {ji,j 2 , • • • ,ji} where ji < j 2 < • • • < ji- If there 
is a loop going through ji,j 2 , • • • ,ji successively, then there will be a term in the product. 




3233 ■ ■ -jljl ^jl32 ^j233 


■^3131 


— ( 1) £23'-/lfljij2®i2i3 ■ ■ ■ ^jljl ~ ( 1) ^k)^jlj2^32j3 ■ ■ ■ 


where a is an element of the symmetric group Si, 


/ 1 2 ••• 1-1 I \ 

\2 3 ■■■ I 1 J ' 


( 37 ) 
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On the right hand side of this equation we have used the signature (parity) function e which 
takes value 1 if the group element is an even permutation and takes value —1 otherwise. 
Obviously e{a) = (—1)^“^ which gives an overall factor —1 in (36). For a generic loop going 
through ji to ji but not necessarily in an ascending order, what is the overall factor? We 
have mentioned that equivalent loops differ only on which vertex is the starting vertex. So 
we can safely set all inequivalent loops to start at ji, then any simple loop can be written as 

Now the problem is how to determine other subscripts of e. To determine the p-th subscript 
we need to know which of these aij has its first subscript being jp and then the second 
subscript would be the p-th subscript of e. We first define another element of the same 
symmetric group, 


and its inverse element 


Through this permutation operator we can find out which aij has its first subscript equal to 
jp, that is jp = So the second script would be = 

j^^^-i^py Therefore the e in (38) can be explicitly calculated 

= erar-i(i)...rar-i(z) = e{TaT-^) = e{a) = (-1)'-^ (41) 

Thus we proved our hrst proposition that = —1. Then (32) simplifies to = 

^{ii 3 k} h=k na(“[®]) come to conclusion in the Proposition 7.1. □ 


1 2 
1 ^2 



*2 

2 


(39) 


(40) 


7.1 A Generic Formula for Superpotential 

From the aforementioned, we observed that the inverse of Ihara zeta function of fully directed 
quivers without self-adjoining loops can be explained by the number of disjoint simple cycles 
in the graph. In this setting, the inverse of Ihara zeta function is therefore a finite polynomial. 
More importantly, the inverse only looks at simple cycles in the quiver, which are simply terms 
in a generic superpotential. Therefore the Ihara zeta function, which is an infinite expansion 
of a finite polynomial, itself simply collects all possible combination of simple cycles of various 
lengths. 

While (25) is the general conclusion, we would like to use this formula to calculate the 
reciprocal of Ihara zeta functions in a special case, where any two inequivalent loops share at 
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least one point, so in (25) there will be no product since there is no disjoint loops according 
to this special setup. Therefore here we would have 

Ck = - [a], (42) 

length k cycles 

That is, Cfc is the negative of the number of all inequivalent simple cycles of length k in the 
corresponding graph. Note that traversing the same vertex sequence but through different 
arrows is considered as traversing different cycles. 


In this perspective, the Ihara zeta function gives the most generic superpotential when 
given information on the quiver alone (of course, in obtaining Calabi-Yau moduli spaces for 
certain theories one needs specific couplings). We now introduce a refinement of the Ihara 
zeta function, i.e, coefficients of the monomials in expansion will be more than just one 
complex variable z but one for each arrow in the quiver. To give some illustrative examples, 
let us consider the dPo and Fq (cones over and Pi x Pi) cases separately. 


For cIPq, we have 9 fields 'hi, .,.,9 and we shall replace the entries in the adjacency matrix 
with these fields: 


/ 0 T 1 + T 2 + T 3 0 

Qp 2 = 0 0 4^4 + ^5 + ^6 

\ T 7 + 4'8 + Tg 0 0 

So we have the inverse of Ihara zeta function as 

C/ = 1 - (Ti + T 2 + T3)(T4 + Ts + 4'6)(^7 + ^'8 + ^ 9 )^^ 


(43) 


= 1 - ('hi'h4'h7 + 'h2'h4'h7 + 4'3'h4'I'7 + Ti'I'5'h7 + 'h2'h5'I'7 + 'h3'h5'h7 

+ 'hi'h6'h7 + 'h2'h6'h7 + 'h3'I'6'I'7 + 'hi'h4'I'8 + 'h2'I'4'h8 + 'h3'h4'h8 

+ T1T5T8 + 4'2^5^8 + ^'3^5^'8 + + ^'2^6^8 + ^3^'6^8 

+ T 1 T 4 T 9 + 4'2^4^'9 + ^'3^4^'9 + ^1^5^9 + ^2^5^9 + ^3^'5^9 

+ TiTgTg + 4'2^6^'9 + 4'3^6^'9)-2:^- 


Here we have all possible cubic superpotential terms that can be constructed from the quiver, 
however the true full superpotential for dPo theory is (e.g., see the catalogue in [38]) 


>Vp2 = 'hi'I'4'1'7 - Ti'h64'8 - 'h2'h4'h9 + 'h2'h5'h8 - 'h3'I'5'h7 + (45) 


Therefore, the expansion of Ihara zeta function in terms of bi-fundamental fields will give 
the most generic superpotential of a specific quiver gauge theory, where the constant 1 in the 
expansion should be removed as well as changing the signs of non-constant in Moreover, 
we can also interpret as the coupling for all the terms in superpotential. The moduli space 
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here would be trivially a point since the F-terms have as many non-trivial constraints as 
there are number of fields. In other words, the Jacobian ideal of the generic superpotential 
is trivial. 


Next, let us consider a Seiberg dual phase for the Fq theory. We have its adjacency matrix 
in the following form: 



(46) 


With this adjacency matrix, we have its Ihara zeta function to be 

6 8 18 10 24 


n=i n=Y n=li n=y n=iy 

where it is obvious that all possible combinations of cubic and quartic terms are generated, 



and upon inverting, we have an infinite polynomial generated from the above finite one. 


7.2 Adjoint Fields and Bi-Directional Arrows 


While most of our Af = 1 gauge theories consist only of directed arrows, there are some 


that have self-adjoint loops (adjoint fields) as well as bi-directional arrows, such as the SPP 
theory and the dual phases of theories. For these, the form of the Ihara zeta function 
are more complicated since we can replace each pair of directed arrows between nodes by a 
single undirected edge and it would be interesting to see the difference between treating these 
edges as directed and treating them as undirected. In the following, we will take SPP case 
as an example to illustrate this difference. 

The SPP quiver [39,40] is shown in Figure 8 from which we can easily read off its adjacency 
matrix, 


/ 2 1 I 



(48) 


A= 1 0 I 
\ 1 I 0 


There are several different ways to view its edges: 

1 All edges are considered to be directed, thus a fully directed graph. We then can still 
use (10) to compute Ihara zeta function. 
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Figure 8: SPP Quiver. All arrows are bi-directional and graphically they can be replaced by 
undirected edges. 

2 All edges are considered to be undirected, then it is an undirected graph 
to utilize an undirected version of formula (2), 

(1 _ ^2^-Tr(Q-I)/2 

Det(I — Az -)- Qz'^) ’ 

to compute the zeta function. 

3 Only the self-adjoint loop is considered as undirected, thence we get a partially directed 
graph. In this case we have to use Eq. (2) itself to compute Ihara zeta function. 


. We then need 

(49) 


In the first perspective, the Ihara zeta function is extremely simple, 

C^{z) = {1 + z){l - 3z). 


(50) 


It can be readily seen that there are only two poles z = —1 and z = 1/3 both being real. For 
the above second case, we have Q = diag{3,1,1} and the zeta function follows immediately. 


C^{z) = {l-zf{l + z){l + z + z^){l-2z + 2z‘^-3z^). (51) 


The poles of this function are plotted in Figure 9 as green and blue points. In the last case, 
Q only counts undirected degrees and P is only responsible for directed edges, so we find 
Q = diag{l, —1, —1} and 


P = 


( ° ^ 
1 0 

VI 1 



Then the resulting zeta function reads 


C\z) = {l-z){l + z){l- 


(52) 


30 


2z-z^ - 2z^). 


(53) 





Im[z] 


1 


0.5 


-0.5 


0.5 


Re[z] 

1 


-0.5 






-1 - 


Figure 9: Poles of lhara zeta function for the SPP quiver. The green and blue points are poles 
when we see all the edges as undirected ones, whereas the red and blue points are poles when we 
take SPP quiver as a partially directed graph. 


The poles in this case are also plotted in Figure 9 as red and blue points. 

From the above example, one can see that the behaviors of poles for lhara zeta function 
are quite different when we treat a graph with self-adjoint loops and bi-directional arrows 
from different perspectives. 


Next we turn to the refined lhara zeta function in which the adjacency matrix A takes 
fields as its entries other than just the number of arrows/edges/loops. For the first case 
viewing all edges as arrows, we have 


A = 



which results in the lhara zeta function 


CsFpi^} = 1 - 20Z - + ^2^'5 + ^' 4 ^ 5 ) 2 ^ - + ^'2^3^6 “ ■ (54) 


Therefore we still have the refined lhara zeta function as the generating function of generic 
superpotentials. However, for the second and third cases, we have to face the problem of 
how to define the Q and P matrices. Since Q is in the power, we would assume it remains to 
be a matrix of numbers, whereas P can be assumed to comprise of either numbers or fields. 
Nonetheless, in either assumption the final lhara zeta function cannot be thought of as a 
superpotential generating function. 
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8 Conclusions and Prospects 


In this paper, we have investigated the Ihara zeta function associated to quiver gauge theories, 
especially under the action of Seiberg duality. Aid by a multitude of explicit examples, we 
have learnt many lessons. 

In the case, the Seiberg dual theories are well-known to be in one-to-one correspondence 
to solutions of the Markov Diophantine equation and we proved that all poles of the Ihara 
zeta functions are on the lines of cubic roots of unity. For x P^ theories, there is also 
an underlying Diophantine equation of a so-called 3-block structure [34,35]. Due to the 
complexity here, we adopted a numerical approach and found that every (even) solution with 
its components being even and less than 800 corresponds to some quiver theory that is a 
Seiberg dual phase. 

Inspired by the graph version of the Riemann Hypothesis, We studied the distribution 
of poles of the Ihara zeta function along the Seiberg duality tree of the various theories, 
including the various del Pezzo and Hirzebruch geometries. The results are shown in figure 
6 where we can see that in all cases, the poles are concentrated in strip areas on or near the 
coordinate axes, which shows an axial symmetry about the horizontal axis. While there is 
yet to be a formulation for the Riemann Hypothesis for irregular digraphs as in our quivers, 
it is still reasonable to speculate that all poles should be constrained within some area other 
than randomly distributed everywhere in the complex plane. 

Nevertheless, we can still use the definitions of strong and weak versions of the graph 
Riemann Hypothesis (cf. Eqs.(5) and (6)) with certain choice of (max) degree q to see how it 
acts on Seiberg duals. Numerical results are summarized in §6.1. It is interesting to observe 
how this RH setup selects certain quivers in the duality tree and hence produce specihc prime 
sequences in the coefficients of the zeta function. 

As a interesting by-product, we find a graph-theoretic interpretation for the reciprocal of 
Ihara zeta function which is summarized in (25) and which has gauge theoretic repercussions. 
As a consequence of this formula, it is immediate to conclude that a refined version of the 
Ihara zeta function is a generating function for the generic superpotential (with couplings in 
powers of z) for the quiver gauge theory. 

One of the immediate challenges is to have a reasonable generalization of the graph 
Riemann Hypothesis for directed quivers. While there is a lack of functional equation in such 
a case, we have seen in our cases that the distribution of the poles is not random and that 
under Seiberg duality, they remain in constrained regions. 


32 


In section 7.2, we have seen that the refined Ihara zeta function cannot be a generating 
function of the SPP quiver. So another interesting work is to find a closed form (may not 
be a Ihara zeta type function) for generic superpotentials of undirected quivers or partially 
directed quivers comprised of directed edges and undirected self-loops. After all, these two 
kinds of quivers both have significant meanings in quiver gauge theories. 
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Appendices 

A Ihara Zeta Function: Determinants and Euler-Product 


In [4] , another zeta function called edge zeta function is constructed to give a determinant 
form of Ihara zeta function for undirected graph. Before giving the full definition of this zeta 
function, we need a few other definitions: 


• The edge matrix W of size 2m x 2m for an undirected graph with m undirected edges 
has entries Wij. The (i,j)-th entry of W, Wij, is a complex variable if the edge e* is 
connected with edge ej with ej ^ e)"^, and the entry is 0 if otherwise. 

• For a closed path C in an undirected graph X written as a sequence of edges C = 
6162 • • • Cs, the edge norm of C is 

Ne{C) = W12W23 ■■ - Wsl- 
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With the above definitions in hand, we can define the edge zeta function as follows 

Ce{W,X)= n il-NEiC))-\ 

[P]GPrime Cycles 

It is clear from this definition that if wtj is set to z G C, we recover the original Ihara zeta 
function such that 

CG{z)=CEiWi,G), 

where Wi is the edge matrix when all non-zero entries set to z. 

Furthermore, we have the following theorem (cf. Chapter 3 of [5]) 

Theorem A.l. 

Ce{W,G) = Det{I-W)-\ 


Proof. By taking logarithm of Euler-product, we have 

-logCi<;(VF,G) = EE -Nsipy, 

[P] i>i ^ 


where Taylor expansion of logarithm is used. As we have 1{P) elements in the prime equivalent 
class [P], hence 


-logCE{W,G) = 


E 

m>l 

f>i 


1 

jm 


p 

l[P]=m 


The sum now is over all prime paths. This is equivalent to sum over all paths of the form 
G = P^ with length jm. It then follows that 


-\ogCE{W,G) = Y,^^Ne{G). 


Since we also gave the following 

V ^,Ne{G) = E -TV(IF™) = -Trlog(/ - IT) = -logDet(/ - IT)-\ 

C ^ ' m>l 


where the first equality can be understood as the trace over IT™' is collecting all cycles of 
various length, then we have Ce{W, G) = Det(/ — W)~^. □ 


Before we give proof to the determinant formula for Ihara zeta function, some matrix 
identities are needed. Proofs of these identities are given in Chapter 3 of [5]. The following 
definitions are needed in the proof: 
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First set J = { ™ V Then the n x 2m start matrix S and the n x 2m terminal 

\Im 0 ) 

matrix T are defined as 


S„p. = 


1, if r; is the starting vertex of edge e, 
0, otherwise, 


0, otherwise. 

With the above definitions, we have the following identities 


T — 

^ve — 


1. SJ = T and TJ = S. 

2. If A is the adjacency matrix of graph G, and Q + In is the undirected degree matrix 
whose j-th diagonal entry specifies the degree of j-th vertex in G. Then we have 
A = ST'^ and Q + In = SS^ = TT^ . 

3. The edge matrix Wi obtained by setting all non-zero entries of VF to 2 has the identity 
Wi + J = T'^S. Here is the usual notation for matrix transpose. 


Now we come to the proof of determinant formula of Ihara zeta function for undirected 
graphs, i.e., Cg(z) = (1 — — Az + This proof is elaborated in 

more details in [4]. First consider a graph with m undirected edges and n vertices, we have 
following matrix equation from previous identities 

/in ^ ^ f '\ _/ In - Az + Qz^ S'z \ / In 0 

V ^2r„ y V 0 - WiZ ) ~ \ 0 hm + JZ ) \ Z hm 

In the above equation, all matrices are of (n -|- 2m) x (n + 2m), and we take determinant on 
both sides of equation to give 

(1 - z2)”Det(/ - Wiz) = Det(4 - Az + Qz^)Det(l 2 m + Jz). 

With the observation that 

implies 

I 

-Iz 


I+ Jz = 
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we obtain Det(I + Jz) = (1 — Therefore it is clear that 

Det(/ — Wiz) = Det(/n — Az + Qz^)r)et{l — 

where the LHS is now (^g{z)~^ from theorem A.l. Since m — n = Tr(Q — -f)/2, we have the 
determinant formula of Ihara zeta function. 

For the details of proof of (2) that gives determinant expression of Ihara zeta function 
for partially directed graphs, see [8]. 


B Solutions to the Markov Equation 


Most of the following proofs can be found in the book [37], we nonetheless rephrase them 
here for the purpose of integrality of our logic flow in Section 4.2. However, in order to 
conform with our underlying physics context, instead of using the original definition of Markov 
equation 

+ v'^ + w'^ = 3uvw , (55) 

we define it here as 

= uvw , (56) 

which is essentially the same. Since in physics we are only interested in positive solutions, 
the positivity of solutions are always presumed in following deductions. 

Lemma B.l. // (a, 6, c) is a solution to eq. (56), {be — a,b,c) is also a solution as well as 
{a, ac — b, c) and (a, b, ab — c). 

Proof. It is readily seen that eq. (56) is symmetric under permutation of variables {u,v,w), 
it is therefore sufficient to prove one of these combinations satisfies eq. (56). 

{be — a)^ + 6^ + = {be — a)he + {a^ + 6^ + — abe) = {be — a)be. 


□ 


Comparing this with the Seiberg duality rules in Section 3.1, one would see that (a, b, e) —)■ 
{be — a, b, e) is exactly of Seiberg duality type apart from a common factor —1. Thus we shall 
call this operation on solutions ^''Seiberg transformation’’'. There are two obvious solutions, 
(3,3,3) and (3,3,6) (of course, (3, 6, 3) and (6, 3, 3) are considered as the same solutions as 
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(3,3,6) due to the symmetry under permutation). These two are dubbed singular solutions 
because there are equal numbers in both solutions. Non-singular solutions are then defined 
as solutions with three distinct numbers. 

Lemma B.2. All solutions to Markov equation except the above two are non-singular. 

Proof. Suppose we have a singular solution (a, b, c), without loss of generality, we shall assume 
a = b, then 2a^ + = a^c => = (c — 2)a^. This indicates \Jc — 2 = c/a is a rational 

number, which means c — 2 has to be a perfect square, say df where d is an integer. Hence 

with c = ad, we have 2a^ + a^d^ = da^ 2 = d{a — d), which means d = 1 or 2. In either 
case we will have a = 3. Then c = ad = 3 or 6 . Thus the two solutions (3, 3, 3) and (3, 3, 6 ). 

Any other solutions would have u, v, w being distinct. □ 


Now for all non-singular solutions, we can write them in a manner such that each has its 
components in ascending order — {u, v, w) with u < v < w. 

Lemma B.3. If we have a non-singular solution (a, b, c) with a < b < c, then ab — c < b and 
be — a > ac — b > c. 


Proof. Consider the quadratic function 

f{x) = — abx + 0 ^ + 1“^ , (57) 

which has two zero points, c and ab — c. Also we can evaluate f{x) at x = b, 

f{b) = — {a — 2)6^ < a^ — 6 ^ < 0 , (58) 

where in the second step we have used the fact that a > 3. Now that f{x) is an upwards 
concave parabola, what immediately follows is that ab — c < b < c. With the same trick we 
can prove ac — b > c and by definition we have be — a > ac — b. □ 

Now assuming a < b < c, we shall call (a, b, c) —)■ {ab — c, a, b) descending transformation 
and the other two ascending transformation. The above lemma tells us that by Seiberg trans¬ 
forming any non-singular solution, we would end up with a “bigger” or “smaller” solution 
depending on which node (or which arrows) your are dualizing on. As for singular solutions, 
if we Seiberg transform (3, 3, 6) we would get (3, 3, 3) (which is smaller) or a non-singular 
solution (which is bigger), and if we Seiberg transform (3, 3, 3) we can only get (3,3, 6 ). Fol¬ 
lowing this logic, if we we have a non-singular solution, then by repeatedly apply descending 
transformation, we will constantly get smaller solutions until there are no smaller ones, that 
is, it will stop at (3,3,3). This also means that we can reverse all transformation steps to 
generate any non-singular solution from (3, 3, 3). Thus we have: 
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Theorem B.4. All positive solutions of Markov equation can he generated from (3, 3, 3) by 
Seiberg transformation. 


The solutions generated from a given solution are called adjacent solutions. For example, 
{3yz — X, y, z) is called the adjacent solution to (x, y, z) with respect to x. 


References 

[1] Y. Ihara, “On discrete subgroups of the two by two projective linear group over p-adic 
fields”. J. Math. Soc. Japan 18 (1966) 219-235. 

[2] M. Horton, “Zeta functions of digraphs”. Linear Algebra Appl. 425 (1) (2007) 130-142. 

[3] H. Mizuno and I. Sato, “Zeta functions of digraphs”. Linear Algebra Appl. 336 (2001) 
181-190. 

[4] H.Bass, “The Ihara-Selberg zeta function of a tree lattice”, Internat.J.Math.,3(1992),TIT- 
TOT. 

[5] Audrey Terras, “Zeta Functions of Graphs: A Stroll through the Garden”, Cambridge 
Studies in Advanced Mathematics (No. 128). 

cf. http: //math.ucsd. edu/~aterras/ 

[6] L-functions in geometry and some applications, Proc. Taniguchi Symp. 1985, “Curvature 
and Topology of Riemannian Manifolds”, Springer Lect. Note in Math. 1201(1986), 266- 
284. 

[7] Fundamental groups and Laplacians, Proc. Taniguchi Symp. “Geometry and Analysis on 
Manifolds”, 1987, Springer Lect. Note in Math. 1339(1988), 248-277. 

[8] Andrei Tarfuleaa, and Robert Perlisb, “An Ihara formula for partially directed graphs. 
Linear Algebra and its Applications, Volume 431, Issues 1-2, 1 July 2009, Pages 73-85. 

[9] R.Murty, “Ramanujan Graphs”, Journal of the Ramanujan Math.Society, 18, No.l. 2007. 

[10] W.C. Winnie Li, “Number Theory with Applications”, Series of University Mathematics, 
Vol. 7, World Scientific, 1996. 

[11] W.C. Winnie Li, “Recent developments in automorphic forms and applications”. Number 
Theory for the Millenium, Volume 2, edited by M.A. Bennett, B.C. Berndt, N. Boston, 
H.G. Diamond, A.J. Hildebrand, and W. Philipp, 2002, pp. 331354, A.K. Peters, Natick, 
Massachusetts 


38 



[12] A. Lubotzky, “Discrete Graphs, Expanding Graphs and Invariant Measures”, Progress 
in Mathematics, 125 Birkhauser, 1994. 

[13] A. Valette, “Graphes de Ramanujan, Asterisque”, 245 (1997), 247276. 

[14] J. McKay, “Graphs, singularities and finite groups”, Proc. Symp. Pure Math. AMS 37: 
183 

J. McKay, “Cartan matrices, finite groups of quaternions, and Kleinian singularities”, 
Proc. AMS, 81, 153 - 154. 

[15] Y. H. He, “Some remarks on the finitude of quiver theories,” [hep-th/9911114]. 

[16] Y. H. He, V. Jejjala and D. Minic, “From Veneziano to Riemann: A String Theory 
Statement of the Riemann Hypothesis,” arXiv:1501.01975 [hep-th]. 

[17] Wai-Kai Chen, Graph Theory and its Engineering Applications, Vol. 5, River Edge, New 
Jersey: World Scientific, 1997. 

[18] Y. H. He, “Graph Zeta Function and Gauge Theories,” JHEP 1103, 064 (2011) 
[arXiv:1102.1304 [math-ph]]. 

[19] N. Seiberg, “Electric - magnetic duality in supersymmetric nonAbelian gauge theories, 
Nucl. Phys. B 435, 129 (1995) [arXiv:hep-th/9411149]. 

[20] K. A. Intriligator, R. G. Leigh and M. J. Strassler, “New examples of duality in chiral 
and nonchiral supersymmetric gauge theories, Nucl. Phys. B 456, 567 (1995) [arXiv:hep- 
th/9506148]. 

[21] P. S. Aspinwall and A. E. Lawrence, “Derived categories and zero-brane stability, JHEP 
0108, 004 (2001) [arXiv:hep-th/0104147]. 

[22] A. Giveon and D. Kutasov, “Brane dynamics and gauge theory, Rev. Mod. Phys. 71, 
983 (1999) [arXiv:hep-th/9802067]. 

[23] A. Hanany and E. Witten, “Type HB superstrings, BPS monopoles, and three- 
dimensional gauge dynamics, Nucl. Phys. B 492, 152 (1997) [arXiv:hep-th/9611230]. 

[24] S. Elitzur, A. Giveon and D. Kutasov, “Branes and N = 1 duality in string theory, Phys. 
Lett. B 400, 269 (1997) [arXiv:hep-th/9702014]. 

[25] B. Feng, A. Hanany and Y. H. He, “D-brane gauge theories from toric singularities and 
toric duality,” Nucl. Phys. B 595, 165 (2001) [hep-th/0003085]. 

[26] B. Feng, A. Hanany, Y. H. He and A. M. Uranga, “Toric duality as Seiberg duality and 
brane diamonds,” JHEP 0112, 035 (2001) [hep-th/0109063]. 


39 



[27] F. Cachazo, S. Katz and C. Vafa, “Geometric transitions and N=1 quiver theories,” 
hep-th/0108120. 

[28] C. E. Beasley and M. R. Plesser, “Toric duality is Seiberg duality,” JHEP 0112, 001 
(2001) [hep-th/0109053]. 

[29] Sergey Fomin, Andrei Zelevinsky “Cluster algebras I: Foundations”, 
arXiv;math/0104151 [math.RT] 

[30] S. Cecotti and C. Vafa, “On classification of N=2 supersymmetric theories,” Commun. 
Math. Phys. 158, 569 (1993) [hep-th/9211097]. 

[31] B. Feng, A. Hanany, Y. H. He and A. Iqbal, “Quiver theories, soliton spectra and Picard- 
Lefschetz transformations,” JHEP 0302, 056 (2003) [hep-th/0206152]. 

[32] S. Pranco, A. Hanany, Y. H. He and P. Kazakopoulos, “Duality walls, duality trees and 
fractional branes,” hep-th/0306092. 

[33] S. Franco, A. Hanany and Y. H. He, “A Trio of dualities: Walls, trees and cascades,” 
Fortsch. Phys. 52, 540 (2004) [hep-th/0312222]. 

[34] S. Benvenuti and A. Hanany, “New results on superconformal quivers,” JHEP 0604, 
032 (2006) [hep-th/0411262]. 

[35] A. Hanany, Y. H. He, C. Sun and S. Sypsas, “Superconformal Block Quivers, Duality 
Trees and Diophantine Equations,” JHEP 1311, 017 (2013) [arXiv:1211.6111 [hep-th]]. 

[36] Li, W.-C. W., “Character sums and abelian Ramanujan graphs”, J. Number Theory, 41 
(1992), 199-217. 

[37] J. W. S. Cassels, “An introduction to Diophantine approximation”, Cambridge Univer¬ 
sity Press (1965). 

[38] J. Davey, A. Hanany and J. Pasukonis, “On the Classification of Brane Tilings,” JHEP 
1001, 078 (2010) [arXiv:0909.2868 [hep-th]]. 

[39] D. R. Morrison and M. R. Plesser, “Non-spherical horizons. I,” Adv. Theor. Math. Phys. 
3 (1999) 1-81, [arXiv;hep-th/9810201]. 

[40] S. Franco, A. Hanany, K. D. Kennaway, D. Vegh and B. Wecht, “Brane dimers and 
quiver gauge theories,” JHEP 01 (2006) 096, [arXiv:hep-th/0504110]. 


40 



